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Abstract 

We present an analytical next-to-leading order QCD calculation of the partonic cross sections for 
the process pp —)■ (jet h)X, for which a specihc hadron is observed inside a fully reconstructed jet. 
In order to obtain the analytical results, we assume the jet to be relatively narrow. We show that 
the results can be cast into a simple and systematic form based on suitable universal jet functions 
for the process. We conhrm the validity of our calculation by comparing to previous results in 
the literature for which the next-to-leading order cross section was treated entirely numerically 
by Monte-Carlo integration techniques. We present phenomenological results for experiments at 
the LHC and at RHIC. These suggest that pp —)■ (jet h)X should enable very sensitive probes of 
fragmentation functions, especially of the one for gluons. 



1 Introduction 


Final states produced at high transverse momentum {pr), such as jets, single hadrons, or prompt 
photons, have long been regarded as sensitive and well-understood probes of short-distance QCD 
phenomena. Recently, a new “hybrid” type of high-p^ jet/hadron observable has been proposed 
and explored theoretically [1-5]. It is defined by an identified specific hadron found inside a fully 
reconstructed jet, giving rise to a same-side hadron-jet momentum correlation. This correlation 
may for example be described in terms of the variable Zh = Pt/Pt^i where pt and are 
the transverse momenta of the hadron and the jet, respectively. The production of identified 
hadrons in jets was first considered for the case of e+e“ annihilation [1-3] and subsequently 
also for pp scattering [4]. Experimental studies have been pioneered in pp —)■ (jet h)X at the 
Tevatron [6]. At the LHC, the ATLAS [7, 8] and CMS [9] experiments have studied pp —)■ (jet h)X, 
and measurements are being carried out by ALICE [10]. Measurements of the cross section (and, 
perhaps, spin asymmetries) should also be possible at RHIC. 

There are several reasons why it is interesting to study the production of hadrons inside 
jets. Perhaps most importantly, the observable provides an alternative window on fragmentation 
functions [4]. The latter, denoted here by D^[z,pi), describe the formation of a hadron h from a 
parent parton c = q,q, g. The variable 2 ; is the fraction of the parton’s momentum transferred to 
the hadron, and fi denotes the factorization scale at which the fragmentation function is probed. 
Usually, fragmentation functions for a hadron h are determined from the processes e’''e“ —)■ hX 
or ep —)■ ehX. The power of these processes lies in the fact that they essentially allow direct 
scans of the fragmentation functions as functions of z. The reason for this is that to lowest order 
(LO) in QCD, it turns out that z is identical to a kinematic (scaling) variable of the process. For 
instance, in —?• hX one has z = 2ph ■ q/q^ to LO, where ph is the momentum of the observed 

hadron and q the momentum of the virtual photon that is produced by the e+e“ annihilation. 
NLO corrections dilute this direct “local” sensitivity only little. A drawback of —?• hX 

or ep —)■ ehX is on the other hand that the gluon fragmentation function can be probed only 
indirectly by evolution or higher order corrections. 

Being universal objects, the same fragmentation functions are also relevant for describing 
hadron production in pp-scattering. So far, one has been using the process pp —>■ hX as a 
further source of information on the D^{z,p) [11-13]. Although this process does probe gluon 
fragmentation, its sensitivity to fragmentation functions is much less clear-cut than in the case 
of e+e“ —)■ hX or ep —)■ ehX. This is because for the single-inclusive process pp —)■ hX the 
fragmentation functions arise in a more complex convolution with the partonic hard-scattering 
functions, which involves an integration over a typically rather wide range of ^ already at LO. 
As a result, information on the D^{z,p) is smeared out and not readily available at a given hxed 
value of 

The process pp —)■ (jet h)X allows to overcome this shortcoming. As it turns out, if one 
writes its cross section differential in the variable Zh introduced above, then to LO the hadron’s 
fragmentation function is to be evaluated ed, z = Zh- This means that by selecting Zh one can 
“dial” the value at which the D^{z,fi) are probed, similarly to what is available in e^e“ —)■ hX 
or ep —)■ ehX. Thanks to the fact that in pp scattering different weights are given to the various 
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fragmentation functions than compared to —)■ hX and ep —)■ ehX, it is clear that pp —)■ 

(jet h)X has the potential to provide complementary new information on the especially 

on gluon fragmentation. Data for pp —)■ (jet h)X should thus become valuable input to global 
QCD analyses of fragmentation functions. At the very least, they should enable novel tests of the 
universality of fragmentation functions. We note that similar opportunities are expected to arise 
when the hadron is produced on the “away-side” of the jet, that is, basically back-to-back with 
the jet [14], although the kinematics is somewhat more elaborate in this case. 

The production of specific hadrons inside jets may also provide new insights into of the struc¬ 
ture of jets and the hadronization mechanism. Varying Zh and/or the hadron species, one can map 
out the abundances of specihc hadrons in jets. Particle identification in jets becomes particularly 
interesting in a nuclear environment in AA scattering, where distributions of hadrons may shed 
further light on the phenomenon of “jet quenching”. Knowledge of fragmentation functions in 
jet production and a good theoretical understanding of the process pp —)■ (jet h)X are also cru¬ 
cial for studies of the Collins effect [15-17], an important probe of spin phenomena in hadronic 
scattering [18]. 

In the present paper, we perform a new next-to-leading order (NLO) calculation of pp —)■ 
(jet h)X. In contrast to the previous calculation [4] which was entirely based on a numerical Monte- 
Carlo integration approach, we will derive analytical results for the relevant partonic cross sections. 
Apart from providing independent NLO predictions in a numerically very efficient way, this offers 
several advantages. In the context of the analytical calculation, one can first of all explicitly 
check that the final-state collinear singularities have the structure required by the universality of 
fragmentation functions, meaning that the same fragmentation functions occur for pp —)■ (jet h)X 
as for usual single-inclusive processes such as pp —)■ hX. We note that to our knowledge this has not 
yet been formally proven beyond NLO. Also, as we shall see, the NLO expressions show logarithmic 
enhancements at high Zh, which recur with increasing power at every order in perturbation theory, 
eventually requiring resummation to all orders. Having explicit analytical results is a prerequisite 
for such a resummation. In Ref. [3], considering the simpler case of e'''e“-annihilation, such 
resummation calculations for large Zh were presented. 

Technically, we will derive our results by assuming the jet to be relatively narrow, an ap¬ 
proximation known as “Narrow Jet Approximation (NJA)”. This technique was used previously 
for NLO calculations of single-inclusive jet production in hadronic scattering, pp —)■ jet X [19- 
23]. The main idea is to start from NLO “inclusive-parton” cross sections for the processes 
ab —?• cX, which are relevant for the cross section for pp —)• hX. They are a priori not suitable 
for computing a jet cross section, which is evident from the fact that the require collinear 
subtraction of final-state collinear singularities, whereas a jet cross section is infrared-safe as far 
as the final state is concerned. Instead, it depends on the algorithm adopted to define the jet and 
thereby on a generic jet (size) parameter TZ. As was shown in Refs. [19-23], at NLO one may 
nonetheless go rather straightforwardly from the single-inclusive parton cross sections to the jet 
ones, for any infrared-safe jet algorithm. The key is to properly account for the fact that at NLO 
two partons can fall into the same jet, so that the jet needs to be constructed from both. In fact, 
within the NJA, one can derive the translation between the dd^^ and the partonic cross sections 
for jet production analytically. We note that the NJA formally corresponds to the limit 77. —)■ 0, 
but turns out to be accurate even at values 77. ~ 0.4 —0.7 relevant for experiment. In the NJA, the 
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structure of the NLO jet cross section is of the form ^log(-R) + i3; corrections to this are of 0{JZ‘^) 
and are neglected. In this paper, we apply the NJA to the case of pp —)■ (jet h)X, using it to derive 
the relevant NLO partonic cross sections. In the course of the explicit NLO calculation, we hnd 
that the partonic cross sections for pp —)■ jet X and pp —)■ (jet h)X may be very compactly formu¬ 
lated in terms of the single-inclusive parton ones convoluted with appropriate perturbative 
“jet functions”. These functions are universal in the sense that they only depend on the type of 
the outgoing partons that fragment and/or produce the jet, but not on the underlying partonic 
hard-scattering function. On the basis of the jet functions, the NLO partonic cross sections for 
pp —)■ jet X and pp —?• (jet h)X take a very simple and systematic form. In fact, it turns out that 
for pp —)■ (jet h)X the jet functions have a “two-tier” form, with a first jet function describing 
the formation of the jet and a second one the fragmentation of a parton inside the jet. We note 
that the concept of jet functions for formulating jet cross sections is not new but was introduced 
in the context of soft-collinear effective theories (SCET) [1-3, 24-26], although applications to 
pp —)■ (jet h)X have to our knowledge not been given. Jet functions in a more general context of 
SCET or QCD resummation have been considered in Refs. [27] and [28], for example. We also 
note that in Ref. [14] the NLO corrections for the case of away-side jet-hadron correlations were 
presented in the context of a Monte-Carlo integration code. 

Our paper is organized as follows. In Sec. 2 we present our NLO calculation. In particular. 
Sec. 2.3 contains our main new result, the formulation of pp —)■ (jet h)X in terms of suitable jet 
functions. Section 3 presents phenomenological results for pp —)■ (jet h)X for LHC and RHIC. 
We hnally conclude our work in Sec. 4. The Appendices collect some technical details of our 
calculations. 


2 Associated jet-plus-hadron production in the NJA 

2.1 Single-inclusive hadron production in hadronic collisions 


Our formalism is best developed by hrst considering the process H 1 H 2 —)■ hX, where a hadron h 
is observed at large transverse momentum pt, but no requirement of a reconstructed hadronic jet 
is made. This is of course a standard reaction, for which the NLO corrections have been known 
for a long time [29, 30]. The factorized cross section at given hadron px and rapidity rj reads 


daH,H,^hX 

dpTdrj 


2pT 


abc 


dXa rHi( \ 

- fa {Xa^PF) 

n. 


Xb 


X 



dZc da^^f,{s,pT,fl,pF,p'F,PR) 
zl vdvdw 


DUzcI^'f) , 


( 1 ) 


with the usual parton distribution functions the fragmentation functions Df, and the hard- 
scattering cross sections for the partonic processes ab —)■ cX', X' denoting an unobserved 
partonic hnal state. Dehning 


1 / = 1 - 


Pt 


z-v 


W = 


Pt 


SV{1-V) 


( 2 ) 


3 





with \/S the hadronic c.m.s. energy, we have 


0 :^“ = W, xr = 


1 - y 


VWjx^ ’ 


= 


1 - 1 / vw 

-+ — 


Xf, 


Xa 


(3) 


The are functions of the partonic c.m.s. energy s = XaXbS, the partonic transverse momentum 
Pt = Pt IZc and the partonic rapidity fj = rj — ^ \og{xalXb). Since only px depends on Zc, the last 
integral in Eq. (1) takes the form of a convolution. The variables v and w in (1) are the partonic 
counterparts of V and W: 


-V 


, Ptg 
n = 1 -^ , w = 


One customarily expresses px and 17 by n and w. 


Px 

sn(l — v) 

(4) 


(5) 


Finally, the various functions in Eq. (1) are tied together by their dependence on the initial- and 
hnal-state factorization scales, px and p'p, respectively, and the renormalization scale pr. 


The partonic hard-scattering cross sections may be evaluated in QCD perturbation theory. We 
write the perturbative expansion to NLO as 


dvdw 


da 


c,(0) 

ah 


dv 


5(1 — w) -|- 


asipn) dd^J^^ 
271 dvdw 




( 6 ) 


where we have used that w = 1 for leading-order (LO) kinematics (since the unobserved partonic 
hnal state X' consists of a single parton), equivalent to 2f)x cosh.{fj)/\/^ = 1. The NLO terms 
have been presented in Refs. [29, 30]. 


2.2 Translation to single-inclusive jet cross section via jet functions 


As shown in Refs. [19-22], one can transform the cross section for single-inclusive hadron pro¬ 
duction to a single-inclusive jet one. References [19-22] explicitly constructed this translation at 
NLO. We may write the jet cross section as 


daH^H2^ietx 

dp^pdrji^^ 


2 ^ 

s 


E 

ab 


Xa 


Xai 


Xb 




vdvdw 


(7) 

where pip and are the jet’s transverse momentum and rapidity, and where TZ denotes a 
parameter specifying the jet algorithm. For the jet cross section we still have 


x„ 


= IV, xf^ = 


1 - 1 / 


1 - VWjx^ ’ 


as in (3), but with V and W now dehned by 

iet 




^1/(1 - V) ■ 


( 8 ) 

( 9 ) 
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Likewise, v and w are as in (4) bnt with pT —)■ p'-p'- Fnrthermore, in analogy with the inclnsive- 
hadron case, p = | log(xa/a;6). We note that the partonic cross sections relevant for 

jet prodnction depend on the algorithm nsed to define the jet. They do not carry any dependence 
on a hnal-state factorization scale. 


In order to go from the inclnsive-parton cross sections to the jet ones dd^^*, the idea is to 
apply proper correction terms to the former. The dd®^ have been integrated over the fnll phase 
space of all hnal-state partons other than c. Therefore, they contain contribntions where a second 
parton in the hnal state is so close to parton c that the two shonld jointly form the jet for a given 
jet dehnition. One can correct for this by snbtracting snch contribntions from dd®j, and adding a 
piece where they actnally do form the jet together. At NLO, where there can be three partons 
c, d, e in the hnal state, one has after snitable snmmation over all possible conhgnrations: 




= [dKb - - ddj^] + [dd, 


c(e)i 


'ab 


da 


d(c) 

ab 


ddtb^] + [dKb - dd'^X’ - dd^^“^] 


e(c) 




de 


dKi + da"^^ + dd„b 


( 10 ) 


Here dd^X is the cross section where parton j prodnces the jet, bnt parton k is so close that it 
shonld be part of the jet, and dd^^ is the cross section when both partons j and k jointly form 
the jet. The decomposition (10) is completely general to NLO. It may be applied for any jet 
algorithm, as long as the algorithm is infrared-safe. As mentioned before, a property of the dd^X 
is that all dependence on the hnal-state factorization scale p'p, which was initially present in the 
ddX mnst cancel. This cancellation comes abont in (10) becanse the dd^ possess hnal-state 
collinear singnlarities that reqnire factorization. This introdnces dependence on p'p in exactly the 
right way as to compensate the /i'^-dependence of the dd^^. 


In the NJA, the correction terms dd^"' and dd^^ may be compnted analytically. At NLO, they 
both receive contribntions from real-emission 2 —)■ 3 diagrams only. For the NJA one assnmes that 
the observed jet is rather collimated. This in essence allows to treat the two ontgoing partons j and 
k as collinear. The relevant calcnlations for the standard conel and (anti-)/Ci [32-34] algorithms 
were carried ont in Refs. [21, 22], while Ref. [23] addressed the case of the algorithm 

proposed in [35, 36]. We note that we always dehne the fonr-momentnm of the jet as the snm of 
fonr-momenta of the partons that form the jet. This so-called recombination scheme” [37] is 
the most popnlar choice nowadays. 


By close inspection of (10), we have fonnd that in the NJA the jet cross section may be cast 
into a form that makes nse of the single-inclusive parton prodnction cross sections ddX 


daH^H2^ietx 

djPpdr]}^^ 


2 ^ 

S 


E 

abc 


— fX {Xa,PF) 

Xa 


- Pf) 

Xb 


X 



dZc dd®b(s,PT,h,hF,/iF,/ifl) J 

zl vdvdw ^ 



p'f 



( 11 ) 


^Here we have in mind primarily the “Seedless Infrared Safe Cone” (SISCone) algorithm introduced in Ref. [31] 
which represents the only cone-based jet definition known to be strictly infrared-safe. However, for single-inclusive 
jet cross sections, the lack of infrared-safety of other cone-type algorithms occurs first at next-to-next-to-leading 
order in perturbation theory and hence is not an issue here. 
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with inclusive jet functions Jq and Jg. We have pT = /Zc, and 2 ;™™ and n, w are now 

as in (3) and (4), respectively. Equation (11) thus states that one can go directly from the cross 
section for single-hadron production to that for jet production by replacing the fragmentation 
functions D’f in (1) by the jet functions Jc- The latter are such that any dependence on 
disappears from the cross section. They depend on the jet algorithm and hence on a jet parameter 
TZ. For the cone and (anti-)fct algorithms IZ is just given by the usual jet size parameter R 
introduced for these algorithms, while for the jet algorithm of [35, 36] we have TZ = Ijy/fizc with 
f3 the “maximization” parameter dehned for this algorithm. In the NJA we generally assume 
7?. -C 1 and neglect 0{TZ‘^) contributions. The jet functions then read explicitly 





fv'p 



J 9 



h-F 



(5(1 -z) 


(hi?) 

27r 


2Cf{1 + z‘^) 


/ log(l 


V 1 


+ Pqq{z) log (A^) 


+ 5{l-z)lf^^ + CF{l-z) 


(5(1 -z) - 


(hi?) 

27r 

Q^(hi?) 

27r 


Pgq{z) log (A^(l - z)^) +Cfz 
4CA(l-z + z^f /log(l-z) 

z V 1 


+ Pgg(z) log (X^)+ 6(1 -Z)lf^° 


27r 


2n4 


<}g\ 


log (A^(l 


-zf)+z{l 



( 12 ) 


where Cp = 4/3, Ca = 3 and Uf is the number of active flavors, and where the LO splitting 
functions Pij{z) as well as the “plus”-distribution are dehned in Appendix A. The dependence 
on the jet algorithm is rehected in the terms 1^*®° and which are just numbers that we also 
collect in Appendix A. 


Equation (11) evidently exhibits a factorized structure in the hnal state for the jet cross section 
in the NJA. Its physical interpretation is essentially that the hard scattering produces a parton 
c that “fragments” into the observed jet via the jet function JT),, the jet carrying the fraction 
Zc of the produced parton’s momentum. At NLO, the factorization is in fact rather trivial. To 
get a clear sense of it, it is instructive to see how one recovers (7),(10) from (11). To this end, 
we combine (6) and (12) and expand to hrst order in the strong coupling. The products of the 
with the LO J(1 — z) terms in J'c just reproduce the single-inclusive parton cross sections at 
pT = pi/*, i.e. the terms daf^^, da^^ in (10). The only other terms surviving in the expansion 
to 0{as) are the products of the LO terms J(1 — w) /dv of (6) with the 0(as) terms in the 
jet functions. These precisely give the remaining contributions — — (plus the other 

combinations) in (10). Because of the convolution in Zc in (11), the (5(1 — t(;)-function in the Born 
cross section actually Exes Zc to the value Zc = 2pl/* cosh(i 7 )/-\/j. Based on our NLO calculation, 
we evidently cannot prove the factorization shown in (11) to beyond this order. We note, however, 
that similar factorization formulas have been derived using Soft Collinear Effective Theory (SOFT) 
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techniques [24-26], for the case of jet observables in e’''e“ annihilation. In particular, functions 
closely related to our inclusive jet functions Jq^g may be found in Ref. [25], where they are termed 
“unmeasured” quark (or gluon) jet functions. We shall return to comparisons with SCET results 
below. 


2.3 Hadrons produced inside jets 


We are now ready to tackle the case that we are really interested in, H1H2 —)■ (jet h)X where the 
hadron is observed inside a reconstructed jet and is part of the jet. Our strategy for performing 
an analytical NLO calculation will be to use the NJA and the same considerations as those that 
gave rise to Eq. (10). Subsequently, we will again phrase our results in a simple and rather general 
way in terms of suitable jet functions. 

The cross section we are interested in is specified by the jet’s transverse momentum and 
rapidity and by the variable 

— Jet ’ 

Pt 

where as in section 2.1 pt refers to the transverse momentum of the produced hadron. As we are 
working in the NJA, we consider collinear fragmentation of the hadron inside the jet. Thus, the 
observed hadron and the jet have the same rapidities, r] = pj®*, since differences in rapidity are 
0(71“^) effects and hence suppressed in the NJA. 


The factorized jet-plus-hadron cross section is written as 




2pj 


5 


E 

a,6,c 


Xa 


Xa,PF) [ —fb^{Xb,PF) 


X 



^^(jetc)^.^^jet^ h) PF, Pf, PR, Zp) j^h 

vdvdwdzp ® 



(14) 


where x™™, x™™, and fj = r^j®* — 1 log(xa/xfe) are as for the single-inclusive jet cross section, and 
where Zp is the partonic analog of Zh- In other words, the are the partonic cross sections for 

producing a final-state jet (subject to a specified jet algorithm), inside of which there is a parton 
c with transverse momentum = Zppip that fragments into the observed hadron. The argument 
of the corresponding fragmentation functions is fixed by pt = zp^ and hence, using (13), is given 
by 2 ; = Zhjzp. Thus the new partonic cross sections are in convolution with the fragmentation 
functions. Note that all other variables V, JE and have the same definitions as in the single- 
inclusive jet case; see Eq. (9). 


At lowest order, there is only one parton forming the jet, and this parton also is the one that 
fragments into the observed hadron, implying Zp = 1. The partonic cross sections hence have the 
perturbative expansions 


dd. 


(jet c) 
ab 


dvdwdzr. 


J-c,(0) / X 7 -(jetc),(l) 


dv 


TT dvdwdzn 


( 15 ) 
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with the same Born terms /dv as in (6). 

In order to derive the NLO partonic cross sections we revisit Eq. (10). Since we now 

“observe” a parton c in the hnal state (the fragmenting one), we mnst not snm over all possible 
hnal states, bnt rather consider only the contribntions that contain parton c: 

dKb - + d&li + da^f,. (16) 


However, for each term we now need to derive its proper dependence on Zp before combining all 
terms. For the terms and dd^}^\ this is trivial since for all of these terms parton c alone 
prodnces the jet and also is the parton that fragments. As a resnlt, all these terms simply acqnire 
a factor d(l — Zp). This becomes different for the pieces dd^f, dd'^l- Following [21-23], in the NJA 
we may write the NLO contribntion to any da^f as 


da 


cd,{l) 

ab 


a 


dvdw 


= —N'ab^K{v,w,e)6{l - w) 


TT 


dZpZp 


^(1 


r^Pc1<{ 


‘"max,algo 


d^U -2e 
-f- ""jet > 


m 


( 17 ) 


jet 


where we have nsed dimensional regnlarization with D = A — 2e space-time dimensions. Eqna- 
tion (17) is derived from the fact that the leading contribntions in the NJA come from a parton K 
splitting into partons c and d “almost” collinearly in the hnal state. We therefore have an nnder- 
lying Born process ab —)■ KX (with some nnobserved recoil hnal state X), whose ZJ-dimensional 
cross section is contained in the “normalization factor” Mab^K, along with some trivial factors. 
The integrand then contains the ZJ-dimensional LO splitting fnnctions P^{z), where the snper- 
script “<” indicates that the splitting fnnction is strictly at 2 ; < 1, that is, withont its 5(1 — z) 
contribntion that is present when c = K. The fnnctions are dehned in Eq. (A.3) in Appendix A. 
The argnment of the splitting fnnction is the fraction of the intermediate particle’s momentnm 
(eqnal to the jet momentnm) transferred in the splitting. In the NJA it therefore coincides with 
our partonic variable Zp. In the second integral in (17) mjet is the invariant mass of the jet. The 
explicit factor in the denominator represents the propagator of the splitting parton K. The 
integral over the jet mass runs between zero and an upper limit minax,aigo; which in the NJA is 
formally taken to be relatively small. As indicated, mmax.aigo depends on the algorithm chosen to 
dehne the jet. We have [21-23] 


m 


2 

max,algo 


min 

< - ^p) 

min(zp, 1 - Zp) 


cone algorithm, 
(anti—)/ct algorithm, 

Jet algorithm. 


( 18 ) 


To make the cross section dd^^^ differential in Zp we now just need to drop the integration over Zp 
in (17). We next expand the resulting expression in £. The integration produces a collinear 
singularity in 1/e. It also contributes a factor (1 — Zp)~^ at large Zp which may be combined 
with the explicit factor (1 — Zp)~^ in (17). In the presence of a diagonal splitting function in the 
integrand we hence arrive at a term {1 — Zp)~^~‘^^, which may be expanded in e to give a further pole 
in 1/e and “plus”-distributions in 1 — Zp. The double poles 1/e^ arising in this way cancel against 
double poles in dd/^f^dd/jf’. The remaining single poles are removed by collinear factorization 
into the fragmentation function for parton c. For non-diagonal splitting functions there are only 









single poles which are directly subtracted by factorization. We note that the original is 

in fact needed both for the cross section with parton c fragmenting and also for the one where 
d fragments. This is reflected in the fact that the Zp-integral in (17) runs from 0 to 1, while for 
dd'^’^^'^ the limit Zp —)■ 0 is never reached as long as Zh > 0. For parton d fragmenting, however, 
we need to use which differs from only by a change of the splitting function. In 

case of a quark splitting into a quark and a gluon, this change is from P^{z) for an observed 
quark to P^g{z) for an observed gluon. Because of P^g{z) = P^{^ — z) one precisely recovers the 
old expression for the inclusive-jet cross sections when all final states are summed over. Likewise, 
if a gluon splits into a qq or gg, the relevant splitting functions P^{z), P^g{z) are by themselves 
symmetric under z -H- 1 — z. 


From this discussion, and combining with Eqs. (15),(16), we obtain to NLO in the NJA: 


dd 


(jet c) 


ab 


dvdwdzr. 


dd' 


ab 


dd' 


c{d) 

ab 


dd' 


c(e) 

ab 


dvdw dvdw dvdw 


(5(1 — Zp) -\- 


dKi 


+ 


dKl 


dvdwdZr, dvdwdz. 


(19) 


Computing and inserting all ingredients of this expression, we find that the cross section may be 
cast into a form that again makes use of the single-inclusive parton production cross sections dd^f^, 
similar to the case of inclusive-jet production in Eq. (11): 


dp^^^drji'^^dzh 



— fa" {Xa, {Xb, fip) 

Xa Xb 

b 


X 


dZc dd^^^{s,PT,fi,fiF,ii'F,i^R) 


vdvdw 







P-'f 


Pf 


, fVR D^, 



( 20 ) 


where and 2 ;™“ are as given in (3), with V and W dehned in terms of jet transverse 

momentum and rapidity. Furthermore, as in (11) we have px = pl^/zc- The (jet algorithm 
dependent) functions /Cc-s-c' are new “semi-inclusive” jet functions that describe the production of 
a fragmenting parton c' inside a jet that results from a parton c produced in the hard scattering. 
For the “transition” g —)■ g we find 


IC \zz X-^n- 

' Pf Pf 


TZpf 


'’.l*/!) = i(l - J)i(l - %) + 


27r 


(5(1 - Zp) \ 2Cf{1 + z^) 


/ log(l 


V 1 


-h Pqq{z) log (A^) -F Cf{1 - z) 


+ ^(1 - z) hcp(l + zD + P„(z,) log (n^)+Cp{l-z,)+T^p(zP , 

( 21 ) 

where Xf^°{zp) is a function that depends on the jet algorithm. Since we will write our new jet 
functions in a more compact form below, we do not present the other functions /Cc^c' here but 
collect them in Appendix B, along with the X^f°{zp). 
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As indicated, JCq^q carries dependence on two (final-state) factorization scales, fi'p and fXp. The 
former is the same as we encountered in the case of single-inclusive jets in Eqs. (11),(12). It was 
originally introduced in the collinear factorization for the single-inclusive parton cross sections, 
but now has to cancel exactly between the and the /Cc-s-c'- As in the case of single-inclusive 
jets, the cancelation of dependence on ^'p is just a result of the fact that we foremost define our 
observable by requiring a jet in the final state. In this sense, n'p is simply an artifact of the way 
we organize the calculation and is not actually present in the final answer. The scale ^p, on the 
other hand, arises because we now also require a hadron in the final state. Technically it arises 
when we subtract collinear singularities from the The logarithms in /i'^ are thus just the 

standard scale logarithms that compensate the evolution of the fragmentation functions at this 
order. We also note that there are two sum rules that connect the inclusive and the semi-inclusive 
jet functions [1, 2]: 

1 

dZp Zp Zp, A, K, /ijj) -|- KLq^g(yZf Zp^ A, K, P-r)] J^qi^z, A, Pr) , 

1 

dzp Zp [Aj(^_^g(^, Zpj A, R, pr) -I- K-g^q(^z^ Zp^ A, r, Pr)] A, pr) . (2^) 

Both are fulfilled by our expressions. Furthermore, dzp}Cq^q{z, Zp] X, k, ^r) reproduces the 
quark splitting contributions to J'g, i.e. the first two lines in Eq. (12). 

We may actually go one step further and decompose the functions /Cc-s-c' into products of jet 
functions that separate the dependence on ^ and Zp. We define two sets of functions: 



jq^q {z,X,hr) = S{l-z) 


27r 


2Cp{l + z^) log (A^) 


+ <5(1 - z)lf^° + Cf{ 1 - z) 

Jq^g {z, A, Pr) = [Pgqiz] log (A2(1 - zf) + Cpz] , 


jg^g {Z,X,fiR) = <5(1 -^) - 


Q;s(Pr) 

27r 


dCAl- 




+ <5(1 - z)/fs° 




2^ 


[Pqg{z) log (A2(1 - zf) + z{l - z)] , 


(23) 
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(where as before A = TZp!^^/fi'p), and 


jq-^q h'-fj) — ^p) "h 


(hi?) 

27r 


2C 


’(l + 4)( ."'’b + PkUp) log (k^) 


+ Cp{l- z,) + Iff-iz,) + i(l - “ 


Jm (%. A'h) = [^!w(%) log ('5^(1 - %)^) + CfZp + . 

;; ,, 'j _ I '^s(hi?) 4C*^(1 — Zp + Zp) ^log(l —Zp)^ ^ ^ 

^p-i-s (^p; hi?) — '^(1 ^p) ^ 1 — z / ^99\^p) ) 


+ X^?°(zp) + 5(1 - zp)If ° 


(^p, K, fiR) = [Pg 3 (^p) log {k\1 - Zpf) + Zp{l - Zp) + Xf^°{zp)\ , (24) 

where again k = TZj^p /fip and the /g'|° are as given in Appendix A for the inclusive-jet case. To 
the order we are considering we then have 


A^c^c'(^)-Z^P) A, K, p/{) ^ ^ Jc^e('^) A, p./^) Je^c'(-^p; hi?) ) 


(25) 


and hence from (20) 

^^iiiii2^(jeth)x 


2p. 




Y1 [ . —fa^Xa^Pp) [ —f^^{Xb,PF) 


^ dxh 


dj^pdrii'^^dzh ^ 


X 


X 


a,h,c 

^ ddlh{s,PT,fi,pp,p'p,pR) 

min Zf, vdvdw 


iet 

rk*^ 


.■ I ^ '^Pt 

/ ^ Jc-^e ( ^ci j ) hi? 


Pf 


yr^ dz 

2_^ I -^3e^c'\Zp, _ p 


P ~A I ~ "^hr ,, \ nh ( // 


'Zh 


Pf 


,/iH /i,M-,h'; • 


(26) 


In other words, in the NJA the production of a jet with an observed hadron factorizes into 
the production cross section for parton c, a jet function describing the formation of a jet 
“consisting” of parton e which has taken the fraction Zc of the parent parton’s momentum, another 
jet function je-s-c' describing a “partonic fragmentation” of parton e to parton d inside the jet, and 
finally a regular fragmentation function D^,. This picture is sketched in Fig. 1. It is interesting to 
see that the structure of the first two lines of Eq. (26) is very similar to that of the inclusive-jet 
cross section (11) when formulated in terms of the jet functions Jc- In fact, if we drop the last 
line and perform the sum over parton-type e, we will exactly arrive at (11), since 


jq^q i.Z, A, hi?) T jq^g (^) A, hi?) dJq{z^ A, hi?) ) 

2npg^q{z,\pR) ^ 3g^g{z,\pR) = Jg{z, X, p r) . (27) 
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Figure 1: Sketch of the production of an observed hadron inside a jet, described in terms of the 
jet functions jc^e aud je^d (see text). 

The last line of (26) thus describes the production of an identihed hadron in the jet. 

We note that at the level of our NLO computation we cannot prove the factorization in (26) 
to all orders. In fact, at 0{as) we can move terms between and je^d- On the other hand, 
it seems very natural that the jet functions that we encountered in the single-inclusive jet case 
should play a role also in this case in the “hrst step” of the formation of the hnal state described 
by the jc-^e- Also, our jet functions je^d are identical to the corresponding functions found in the 
SOFT study [3] of hadrons in jets produced in e+e“-collisions, except for endpoint contributions 
cx 5(1 — Zp) that are necessarily different in the SCET formalism due to the presence of a soft 
function. 

We hnally note that the cross section (26) may also be expressed in terms of the hadron 
kinematics, using the relation 

dpTdrjdzh Zh dpf^^drji^^dzh V 

3 Phenomenological results 



We now present some phenomenological results for associated jet-plus-hadron production. First, 
we compare our analytical calculation in the NJA with the one of [4], where the NLO cross 
section was obtained numerically by Monte-Carlo integration techniques. As in that paper, we 
consider the case of charged hadrons produced in pp collisions at the LHC with center-of-mass 
energy \/~S = 8 TeV. We dehne the jet by the anti-kt algorithm with jet parameter R = 0.4. The 
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Figure 2: Comparison of our results in the NJA to the ones of [4] for LHC kinematics. 


renormalization and initial-state factorization scales are set equal to the transverse momentum of 
the jet, p.R = p-F = Pt^i while the hnal-state factorization scale is chosen as p'p = . The latter 

choice serves to sum logarithms of R to all orders [38, 39], although this only becomes necessary 
for jet sizes much smaller than R = 0.4. As in [4] we use the CTEQ6.6M parton distributions [40] 
and the “de Florian-Sassot-Stratmann” (DSS07) fragmentation functions of Ref. [11]. Our results 
refer to (summed) charged hadrons, i.e. h = hC + h~. 


In Fig. 2 we show the ratio of the cross section in the NJA with that obtained numerically in 
Ref. [4]. The ratio is shown as function of Zh, where the cross sections have been integrated over 
\t]\ < 1 and 30 GeV < px < 200 GeV in hadron rapidity and transverse momentum. As one can 
see, the agreement of the two NLO calculations is very good. The deviations are smaller than 
3% everywhere, which demonstrates the good accuracy of the NJA. We note that in [4] a closely 
related variable Zh is considered, which is dehned as 


—* —♦ j6t 

^ Pt - Pt 

IpVV 


( 29 ) 


This dehnition differs from (13) only by 0{TZ‘^) corrections, which are anyway neglected in the 
NJA. In the limit Zh —)■ 1, the two dehnitions become equivalent. This explains why the ratio in 
Fig. 2 is even closer to unity for larger values of Zh- The excellent accuracy of the NJA observed 
in the hgure is consistent with similar comparisons for the case of single-inclusive jet production 
in the NJA [21, 22]. 
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Figure 3: LO (dashed) and NLO (solid) cross sections for pp —)■ (jet 7r)X for ALICE conditions, 
as functions of pion pt- The bands show the scale dependence of the cross section for variations 
of the scale between pi^^/2 (upper end of bands) and (lower end of bands). The factorization 
and renormalization scales have all been set equal and varied simultaneously. 


Next, we show some results for the kinematics relevant for the ongoing studies in ALICE [10]. 
We consider pp collisions at \/S = 7 TeV and fragmentation into charged pions (tt = 7r++7r“). For 
the rapidity interval we choose |p| < 0.5 and we restrict the jet transverse momentum to 15 GeV 
< pl^* < 20 GeV. As before, the jet is defined by the anti-fc^ algorithm with R = 0.4. We now use 
more modern sets for the parton distributions, GTIO [41], and fragmentation functions, DSS14 [12]. 
All scales are set equal to the transverse momentum of the jet, = p,p = p'p = p^^ = p. Figure 3 
shows the LO (dashed) and NLO (solid) cross sections for associated jet-plus-pion production 
differential in the transverse momentum of the pion. Note that the variable Zh is determined as 
Pt/p^t hence is varied upon integration over pl^*. The bands show the changes of the cross 
sections when the scales are varied in the range jfp /2 < p < 2p!p'. As one can see, the scale 
dependence of the cross sections improves somewhat when going from LO to NLO, although not 
as much as one would have hoped. This feature was also observed for single-inclusive hadron 
production in hadronic scattering [30]. 

For the same kinematical setup we also show the cross section differential in Zh, see Fig. 4. A 
fixed value of Zh implies that the hadron’s transverse momentum varies as we integrate over pp^. As 
discussed in the Introduction, this arguably is the most interesting distribution for pp —)■ (jet 7i)X 
since it allows direct scans of the fragmentation functions. Apart from the scale variation, we also 
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Figure 4: Same as in Fig. 3, but as function of Zh- As before the solid bands show the scale 
uncertainty. The hatched band displays the uncertainty of the cross section related to the frag¬ 
mentation functions. This hand is only reliable up to Zh = 0.65 and extrapolated beyond (see 
text). 

show in the hgure the uncertainty related to the fragmentation functions, which we compute using 
the Hessian error sets provided in the DSS14 set [12], Note that the resulting uncertainty band 
is reliable only up to Zh ~ 0.65, beyond which there are presently hardly any hadron production 
data available for annihilation or ep scattering. We hence stop the main uncertainty band 
there and only sketch its possible extrapolation to higher Zh- It is clear from the hgure that 
precise measurements of the cross section as a function of Zh have the potential to provide new 
information on fragmentation functions that is complementary to - and in some respects better 
than - that available from e’''e“ annihilation. 

An interesting question is of course which of the fragmentation functions are primarily probed 
when the cross section for pp —)■ (jet 7r)X is studied as a function of Zh- Depending on kinematics, 
different initial-states may dominate the contributions to the cross section, resulting also in dif¬ 
ferent weights with which the fragmentation functions for the various parton species enter. Given 
how little information on gluon fragmentation is available from —)■ hX and ep —)■ ehX, it 

is especially interesting to see how strongly the cross section for pp —)■ (jet7r)X depends on Dg. 
It is known that for LHC energies, channels with gluonic initial states (especially gg) typically 
make important contributions to cross sections. In order to explore whether this allows probes of 

at the LHC, we investigate in Fig. 5 the relative contributions of quark/antiquark (summed 
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Zh 


Figure 5: Normalized quark (solid) and gluon (dashed) contributions to the cross section differ¬ 
ential in Zh for the kinematic conditions chosen for Fig. 4. We show results for DSS07 [11] and 
DSS14 [12] fragmentation functions. 


over all flavors) and gluon fragmentation to the cross section for pp —)• (jet 7r)X at ALICE (as 
shown in the previous Fig. 4). We normalize the contributions to the full cross section, so that 
the quark and gluon contributions add up to unity. We use both the DSS07 and DSS14 sets. 
As one can see, for Zh ^ 0.5 the two sets give similar results and show that the cross section is 
strongly dominated by gluon fragmentation here. This is already interesting, since it implies that 
in this regime clean probes of should be possible that should be much more sensitive than 
e’''e“ annihilation. Beyond Zh = 0.5, the two sets of fragmentation functions show very different 
behavior. For DSS07, gluon fragmentation continues to dominate all the way up to Zh ^ 0.9, 
whereas for DSS14 the quarks take over at Zh ~ 0.7. We stress again that the uncertainties of the 
fragmentation functions become very large at such values of z, as we saw in the previous hgure, 
and are in fact hard to quantify reliably. It is evident that information from pp —)■ (jet n)X in 
this regime will be most valuable, regardless of whether quark or gluon fragmentation dominates. 
Detailed measurements for various bins in transverse momentum and rapidity will likely help in 
disentangling fragmenting quarks and gluons. 

As mentioned in the Introduction, measurements of charged hadrons produced in jets are 
already available from ATLAS [7] and CMS [9]. ATLAS has published measurements at y/S = 
7 TeV [7] and presented preliminary data [8] also at y/S = 2.76 TeV. The two analyses each 
use a slightly different dehnition of Zh which however both coincide with our Zh in the NJA 
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Figure 6: NLO cross section for pp —)■ {jeth)X as function of Zh at y/S = 7 TeV, compared to the 
ATLAS data [7] for charged hadron production in the leading jet. The cross section is normalized 
to the total jet rate. In the region outside the validity of the DSS07 set the theory curves are 
extrapolated and plotted as dotted lines. 


limit. Figures 6 and 7 present comparisons of our NLO calculations to the ATLAS data for 
the two energies. We have now gone back to the DSS07 set, since unspecihed charged-hadron 
fragmentation functions are not available in the more recent DSS14 set. As one can see, there is 
overall a very good agreement. Note that this agreement extends even down to values of Zh < 0.05, 
well outside the region of validity of the DSS sets. The hgures clearly demonstrate the potential 
of the data to further pin down the charged-hadron fragmentation functions. 

The CMS analysis [9] starts from a dijet sample and then studies charged-hadron production 
inside either the leading jet (which is required to have p!^^ > 100 GeV) or the subleading jet (with 
p!p > 40 GeV). As such, these conditions are different from the single-inclusive jet situation we 
consider in this paper, and strictly speaking we cannot compare to the CMS data. On the other 
hand, it turns out that the CMS data for hadron production in the leading and the subleading 
jet are in remarkable agreement for Zh > 0.05, when one normalizes each of them individually to 
the corresponding total (leading or subleading) jet event rate. This finding clearly indicates that 
fragmentation inside jets is really independent of the underlying event topology and happens in 
the same way in any jet. Therefore, the overall reservation notwithstanding, we show in Fig. 8 the 
comparison of the normalized one-jet rate differential in \og{l/zh) to the CMS data for hadron 
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Figure 7: Same as Fig. 6 but at y/S = 2.76 TeV, compared to the preliminary ATLAS data [8]. 


production in the leading jet. We show the theoretical curve down to Zh ~ 0.1. As one can see, 
the agreement with the data is very good in this regime. We have found that quark and gluon 
fragmentation contribute roughly in equal parts to the cross section. 

We finally note that measurements of pp —>• (jet 7r)X should readily be feasible at RHIC, 
especially in the STAR experiment where both inclusive jet [42] and pion cross sections [43, 44] 
have been measured. Figure 9 shows our NLO predictions as functions of Zh for pp collisions 
at y/S = 200 GeV and y/S = 510 GeV. For the former, we have integrated the jet transverse 
momentum over 5 GeV < < 40 GeV, while for y/S = 510 GeV we have used 10 GeV < < 

80 GeV. In both cases we integrate over |p| < 1. The jet is dehned by the anii-kt algorithm 
with R = 0.6. As before we use GTIO and DSS14 and set all scales equal to the jet transverse 
momentum. 


4 Conclusions and outlook 


We have considered the process pp —)■ (jet h)X, for which a specific hadron is observed inside 
a fully reconstructed jet. Using the approximation of relatively narrow jets, we have performed 
an analytical next-to-leading order calculation of the partonic cross sections for this process. We 
have found that the NLO partonic cross sections may be systematically formulated in terms of 
simple jet functions for the process. These functions are universal; that is, they only depend 
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Figure 8: NLO cross section for pp —)■ {jeth)X differential in ^ = \og{l/zh) at y/S = 2.76 TeV, 
compared to the CMS data [9] for hadron production in the leading jet. The cross section is 
normalized to the total jet rate. 


on the types of partons producing the jet and fragmenting into the observed hadron. We note 
that in the process of computing the jet functions we needed to perform subtractions of the 
final-state collinear singularities. These take the same form as the corresponding subtractions 
in single-inclusive hadron production (without a reconstructed jet). This demonstrates that the 
fragmentation functions are universal to NLO in the sense that the same functions appear in 
pp —)■ (jet h)X as in pp —)■ hX. Essentially, all effects of the fact that a jet is reconstructed along 
with the hadron factorize into a perturbatively computable factor, the jet function. The factorized 
structure in terms of jet functions we hnd at NLO suggests that this statement is true to all orders. 
Our Ending is in line with the result of [3]. 

Our numerical results are in very good agreement with those obtained by Monte-Carlo inte¬ 
gration techniques in [4]. We have presented phenomenological results for the NLO cross section 
for the kinematics relevant for forthcoming measurements at ALICE and for previous ones by 
ATLAS and CMS. These results show that pp —)• (jet h)X should enable very sensitive probes of 
fragmentation functions. In particular, the cross section difierential in Zh probes the fragmentation 
functions almost “locally” at the momentum fraction Zh. The combination of fragmentation func¬ 
tions that is probed depends on the mix of initial-state of parton distributions and hard-scattering 
functions that dominates. We find that, in contrast to the standard process e+e“ —)■ hX that 
is customarily used for extractions of fragmentation functions, the process pp —)■ (jet h)X should 
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Figure 9: NLO cross section for pp —)■ {jeth)X differential in Zh, for STAR kinematics with 
y/S = 200 GeV (dashed) and y/S = 510 GeV (solid). 

offer detailed insights into gluon fragmentation. Also, information at very large Zh might become 
accessible, although here it may become necessary to perform resummations of large logarithmic 
terms in the jet functions. We note that at high Zh typical particle multiplicities in the jet a 
very low, so that power corrections and non-perturbative phenomena will become important here 
as well. As has been discussed in Ref. [45], hadronization corrections to inclusive-jet production 
may exhibit a scaling with 1/R, making them especially relevant in the case of rather narrow jets. 
Although these corrections are at the same time suppressed by an inverse power of transverse 
momentum, it will be an interesting and important task to investigate their structure in case of 
the hadron-plus-jet observable where two separate transverse momenta are present. 

There are various other possible extensions of our work that we hope to address in the future. 
As is well known, hadron production in jets has important applications in studies of spin phe¬ 
nomena in QCD in terms of the Collins effect [18], where the azimuthal distribution of a hadron 
around the jet axis is considered. Studies of the effect in pp scattering [15-17] will require a 
detailed theoretical understanding of the process, to which we hope we have contributed in this 
paper by computing the NLO corrections for the denominator of the spin asymmetry. We expect 
that our method based on jet functions is also applicable to the spin-dependent case. Finally, we 
mention that also photon fragmentation in jets could be interesting as a means to constrain the 
poorly known photon fragmentation functions (see [46] for related work on e’''e“ annihilation and 
ep scattering). 
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A Details for jet functions in the single-inclusive case 

In our results (12) for the single-inclusive jet functions we have the standard LO splitting functions 




1 -l- z'^ 3 . , 

+ -5(1-^) 


^gg(^) — Cp 


1 + (1-^)' 


p - or (1 - ^ + I /^o... X 

Pgg{z) - 2 Ca +2 <^(1 . 

Pqgi^) = ^ ( 2 ^^ + (1 - zf) ) 


(A.l) 


with (3o = — 1^/- The “plus”-distribution is dehned as usual by 


dz f{z)[g{z)]+ = / dz{f{z) - f{l))g{z) 


(A.2) 


Dropping the ^-function contributions and ignoring the “plus”-distributions, we obtain the split¬ 
ting functions P^{z) at z < 1. For our calculations, we actually need these functions computed in 
dimensional regularization in D = 4 — 2e dimensions, where they are denoted as P^{z). We have 






(A.3) 


with 


= -Cp{l - z ), P^fiz) = -Cpz , 
Pig\z) = -Z{ 1 - Z), P^l\z) = 0 . 


(A.4) 


We note in passing that the pieces in (12) that are independent of the jet algorithm may be 
constructed following a simple rule: Each of the jet functions JT), contains the combination 




F.jwiogp"(i-r)-U’w 


(A.5) 
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up to regularization by distributions at z = 1. 


The algorithm-dependent terms and in Eq. (12) may be determined from the calcu¬ 
lations presented in [22, 23]. For cone algorithms we have 


rcone 


= C'F|-^ + y-31og2]) , 

^ , 137 TT^ 11 \ Uf f23 4, 


while for the (anti-)fcj algorithms 




rfe ^ I 67 27r2\ 23 


Finally, for the algorithm; 


^JErr „ / _ vr^ 3 


= Cf ( -5+ y - -log2 ) , 

^ / 45 y 11 Uf /23 2, ^ 

h - C’d-j + y-yiogsl + ylj^ + poga 


Interestingly, we find 


jJej, _ J- / Tcone < Tkt'\ 

~ o V-'i -'i / 


(A.6) 


(A.7) 


(A.8) 

(A.9) 


B Jet Functions for the semi-inclusive case 


In addition to in Eq. (21) we have: 


ICffJz,Zp,X,K,fiR) = 


- 5(1 - Zp)[Pgq{z) log (A^(l - zf) + Cpz 


2e 


<^(1 - ^)\Pgqi^p) log ('^^(1 - ^pf) + CpZp + Xf^^Zp) 


(B.l) 


lCf^Az,Zp,X,K,iiR) = 5(1 - z)5(l - Zp) + 


Qs(/^J?) 

2e 


N - Z + Z^)^ /log(l-z) > n / M /\2\ 

--( ) +^..(^)log(A') 


+ 3(1 - +P„(jJlog(K2)+i;j!"(z,)|| , (B,2) 
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O^sifJ^R) 

2tt 


- <^(1 - Zp)^Pqg{z) log (A2(1 - zf) + z{l -z)^ 


5(1 - z)\Pqg{zp) log {k\ 1 - Zpf) + Zp{l - zp)+lfg^°{zp) 


(B.3) 


where as before A = /^'p, k = Pp^p /^p, and where the algorithm dependent terms are 

2Pc'c{z) log (j^) 0(1/2 — z) cone algorithm, 

2Pc'ciz)\ogz (anti—)A;t algorithm, (B-4) 

Pc'ciz) [log(z) + log (0(1/2 - z)] Jet algorithm. 




We note that a closely related resnlt for the cone algorithm was obtained in [3]. Again, similar 
to (A.9), we have 

d?(.-) = . (B.5) 
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